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SGS  Flux  Vectors  of  Mass  and  Energy 


LIQIU  WANG 

Department  of  Mechanical  Engineering 
The  University  of  Hong  Kong,  Hong  Kong 


Abstract.  Both  necessary  and  sufficient  conditions  are  derived  in  a system- 
atic, rigorous  way  for  a subgrid-scale  (SGS)  flux  vector  model  to  preserve 
the  frame-indifference  of  the  vector  and  to  satisfy  both  the  principle  of  ma- 
terial frame  indifference  (PMFI)  and  the  second  law  of  thermodynamics. 
This  leads  to  the  results  either  confirming  the  previous  intuitive  arguments 
or  offering  new  insights  into  turbulence  modelling,  and  is  of  significance  in 
clarifying  some  controversies  in  the  literature,  examining  how  well  existing 
models  preserve  the  physics,  and  developing  new  models. 


1.  Introduction 

SGS  stresses  and  fluxes  of  mass  and  energy  are  believed  to  be  quantities 
determined  by  filtered  large-scale  velocity  and  mass  fraction/temperature 
fields  in  the  large  eddy  simulation  (LES).  Based  on  this  fundamental,  in- 
trinsic belief,  various  approaches  have  been  proposed  to  relate  SGS  stresses 
and  fluxes  to  the  filtered  large-scale  fields,  so-called  SGS  turbulence  mod- 
elling. The  readers  are  referred  to  Ciofalo  (1994),  Mason  (1994),  Lesieur 
&;  Metais  (1996)  and  Sagaut  (2001)  for  some  recent  excellent  reviews  and 
discussions  of  this  important  topic.  While  some  LES  results  based  on  some 
commonly  used  models  seem  encouraging,  they  fail  to  meet  either  one  or 
both  of  two  natural  fundamental  requirements  for  turbulence  models:  pre- 
serving the  fundamental  properties  of  the  quantities  being  modeled  and 
satisfying  some  classical  principles. 

The  modelling  of  the  SGS  flux  vectors  of  mass  and  energy  consists  of 
replacing  them  by  constitutive  equations  expressing  them  as  functions  of 
filtered  large-scale  fields  of  velocity  and  mass  fraction/temperature.  While 
such  constitutive  equations  may  take  different  forms  such  as  algebraic  and 
differential,  it  appears  to  be  a basic  requirement  to  preserve  the  properties 
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which  the  flux  vectors  hold  by  their  definition.  Such  a property  is  the  frame- 
indifference  (Fureby  & Tabor  1997,  Ghosal  1999,  Wang  2001).  It  follows 
from  the  definition  of  the  SGS  flux  vectors  and  states  that  they  remain  the 
same  directed  line  element  under  a change  of  frame.  The  issue  concerned 
with  whether  a model  guarantees  this  property  is  referred  as  the  invariance 
in  the  literature. 

While  the  first  requirement  focuses  on  the  properties  of  the  SGS  flux 
vectors  themselves,  the  second  requirement  emphases  on  their  function  rela- 
tion with  the  filtered  large-scale  fields.  Such  function  relations  are  required 
to  satisfy  some  classical  principles  including  the  PMFI  and  the  second  law 
of  thermodynamics  (Fureby  & Tabor  1997,  Ghosal  1999,  Wang  1997,  1999, 
2001).  The  PMFI  requires  that  the  function  relation  is  the  same  for  every 
observer,  i.e.  in  every  frame  of  reference.  The  second  law  of  thermodynam- 
ics, on  the  other  hand,  states  that  the  flux  is  always  from  high  concentration 
to  low  concentration.  Note  that  the  realizability  for  the  Reynolds  and  SGS 
stresses  also  comes  from  the  second  law  of  thermodynamics  (Wang  1999, 
2001). 

The  motivation  for  the  present  work  comes  from  the  desire  to  derive 
both  necessary  and  sufficient  condition  in  a systematic,  rigorous  way  for  a 
SGS  flux  model  to  preserve  the  frame-indifference  of  SGS  flux  vectors  and 
to  satisfy  both  the  PMFI  and  the  second  law  of  thermodynamics.  Unlike 
the  works  in  the  literature,  no  intuitive  assumption  is  introduced  in  the 
derivation;  the  independent  variables  are  chosen  properly;  the  PMFI  and 
the  frame  indifference  of  SGS  flux  vectors  are  clearly  distinguished.  This 
leads  to  some  conclusive  results.  Among  them,  some  confirm  the  previ- 
ous intuitive  arguments,  and  others  form  new  insights  to  SGS  turbulence 
modelling. 

2.  Principle  of  Material  Frame- Indifference  and  Second  Law  of 
Thermodynamics 

Consider  a class  of  constitutive  relations  which  relate  the  passive  SGS  flux 
vector  q of  mass  or  energy  to  its  arguments  0, 0.P.,  S/0,  v,  L,  i.e., 

q = f(0,aP.,V0,v,L).  (1) 

Here  f is  a vector-valued  function.  0 is  the  concentration  of  a property. 
It  is  the  mass  fraction  of  a species  for  the  case  of  SGS  flux  of  mass,  and 
the  temperature  when  q is  the  SGS  flux  of  energy.  O.P.  denotes  the  other 
scalar-valued  thermophysical  parameters  which  are  independent  of  v and 
L and  are  typically  the  local  thermodynamic  state  variables,  s/0  is  the 
gradient  of  0.  v is  the  filtered  velocity  vector.  L is  the  velocity  gradient 
tensor  of  v,  a second  order  tensor-valued  variable. 
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In  sharp  contrast  with  that  in  the  literature,  we  choose  L as  an  indepen- 
dent variable  instead  of  its  symmetric  part  D (the  velocity  strain  tensor) 
and  skew  part  W (the  vorticity  tensor)  because  D and  W can  not  be  re- 
garded as  independent.  We  do  not  include  k,  l or  e as  the  independent 
variables.  The  exclusion  of  the  explicit  dependence  of  q on  time  t and  po- 
sition vector  r comes  from  the  fact  that  they  affect  q through  9 , O.P.,  yd, 
v and  L. 

The  relation  (1)  satisfies  both  principle  of  determinism  and  principle 
of  local  action  since  we  assume  that  q at  a point  and  a time  instant  is  a 
function  of  its  arguments  at  that  point  and  that  instant. 

The  principle  of  frame- indifference  requires  that  f is  the  same  for  every 
observer,  i.e. 

q*  = f(r,O.P.\MT,v*,L*)  (2) 

in  which  superscript  * represents  the  quantities  observed  by  another  ob- 
server *. 

The  second  law  of  thermodynamics  states  that  q is  always  from  high 
concentration  to  low  concentration.  This  requires  that:  (1)  f changes  its 
sign  if  y 9 changes  the  sign,  i.e., 

f (9, 0.P.,  - V 9,  v,  L)  = -f (9, 0.P.,  y0,  v,  L),  (3) 

and  (2)  the  projection  of  the  f on  y 9 is  negative  semi-definite,  i.e., 

f • V0  < 0.  (4) 


2.1.  NECESSARY  CONDITIONS  FOR  REQUIREMENTS  (2)  AND  (3) 

2.1.1.  q — v relation 
Theorem  1.  f is  independent  of\. 

Proof  From  the  principle  of  observer  transformations  (Geankoplis  1983, 
Truesdell  1991), 

9*  = 9,  ( O.P.y  = O.P. , q*  = Q(t)q, 

(V0*r  = Q(*M,  L*  - Q(f)LQT(f)  + Q(f)QT(f), 
r*  = Q(t)r  + c(t),  v*  = yf  = Q(*)r  + QWV  + *(*)> 

(5) 

where  Q is  an  arbitrary  rotation  tensor,  r a position  vector  of  material 
point,  c (t)  an  arbitrary  vector-valued  function  of  time  t,  and  a dot  over  a 
letter  indicates  a time  derivative.  In  (5),  we  have  used  the  frame  indifference 
of  q. 
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By  making  use  of  (1)  and  (5),  (2)  yields  (suppressing  t) 

i(6, O.P.,  Qv0,  Qr  + Qv  + c,  QLQt  + QQr)  = Qf (0, 0.P. , y0,  v,  L), 
VQ  and  c. 

(6) 

Since  (6)  holds  for  all  Q,  it  must  be  true  for  Q = 1.  Take  Q = 1,  then 
Q = 0.  Equation  (6)  reduces  to 

f (0, 0.P.,  v$,  v + c,  L)  = f (0, O.P.,  v#>  v,  L)  Vc.  (7) 

This  implies  that  f is  independent  of  velocity  v. 

By  applying  Theorem  1,(1)  and  (6)  reduce  to 

q = f(0,O.P,V0,L),  (8) 

f (0, 0.P.,  Qy0,  QLQr  + QQt)  = Qf (0, 0.P. , y 0,  L)  VQ.  (9) 

2.1.2.  q — L relation 

Note  that  L can  be  uniquely  decomposed  into  a symmetric  tensor  D (ve- 
locity strain  tensor)  and  a skew  tensor  W (vorticity  tensor).  Expression  (9) 
may,  then,  be  rewritten  as 

f(0, O.P.,  Qy0,  QDQT+QWQr+QQT)  = Qf(0, 0.P.,  y0,  L)  VQ.  (10) 

Theorem  2.  For  rotation  tensor  Q (t)  = exp[Cl{t  — r)]Q,  we  can,  at  any 

instant  r,  pick  Q(r)  and  Q(t)Qt(t)  to  be  arbitrary,  independent  rotation 
and  skew  tensors,  respectively.  Here  Q is  any  time-independent  rotation 
tensor,  and  fi  any  time-independent  skew  tensor. 

Proof  As  A is  a time-independent  skew  tensor,  exp[£l(t  — r)]  is  thus 
a rotation  tensor  for  any  fixed  time  r and  all  time  t.  Since  both  Q and 
exp[Cl(t  — t)]  are  rotation  tensors,  Q (t)  = exp[Ct(t  — t)] Q is  also  a rotation 
tensor  for  all  time  t.  Also, 

Q(r)  = Q,  (11) 

Q(r)Qr(r)  = 6Q(t)Qt(t)  = 0.  (12) 

They  are  clearly  independent  rotation  and  skew  tensors  if  Q and  are  any 
time-independent  rotation  and  skew  tensors,  respectively. 

Theorem  3.  L affects  q only  through  velocity  strain  tensor  D. 

Proof  To  prove  this,  choose  Q (t)  defined  in  Theorem  2 as  the  rotation 
tensor  in  (10)  while  for  any  instant  r,  -QWQr  |T  is  used  as  the  skew 
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tensor  fi,  i.e.  ft  = -QWQ7  |r  (Such  aOdoa  skew  tensor  since  fl1  = 
— QW7Qt  |r=  QWQr  | T=  — ft).  Then  at  time  t = r,  (10)  yields 

f (0, 0.P.,  Qy0,  QDQt)  = Qf  (0, 0.P.,  y0,  L)  VQ.  (13) 

As  this  is  true  for  all  rotation  tensor  Q,  it  must  hold  for  Q = 1.  Let  Q = 1, 
(13)  yields 


f(0,O.P.,V0,L)=f(0,O.P.,y0,D), 

(14) 

q = f(0, 0.P.,  v0,  D). 

(15) 

2.1.3.  q-y0  relation 

Expression  (15)  and  the  principle  of  frame-indifference  together  yield 

q*  =f(0*,(O.P.)*,(v0*)*,D*)  (16) 

By  making  use  of  (5),  (15)  and  D*  = QDQr  (Truesdell  1991),  (16)  leads 
to 

f (0, 0.P. , Qy0,  QDQt)  = Qf (0, 0.P.,  V0,  D)  VQ.  (17) 
Also  the  second  law  of  thermodynamics  [Eq.(3)]  requires  that 

f (0, 0.P.,  - V0,  D)  = -f  (0, 0.P.,  V0,  D).  (18) 

Since  the  velocity  strain  tensor  D is  a real,  symmetric  tensor,  it  has 
three  real  eigenvalues.  The  three  eigenvalues  can  be  distinct,  identical,  or 
two  of  them  can  be  identical.  In  the  present  work,  we  focus  on  the  case  that 
the  three  eigenvalues  are  distinct.  Similar  results  may  be  obtained  for  the 
other  two  cases. 

Theorem  4.  \/6,  D\/9,D2\/9  are  linearly  independent  if  three  eigenvalues 
of  D are  distinct. 

Proof  Let  and  ffc  ( k = 1,2,3)  to  be  the  eigenvalues  and  eigenvectors 
of  D.  D may  be  represented,  in  its  spectral  form,  as 

3 

D = ® ffc-  (19) 

k=  1 

The  linear  independence  of  f/,-(fc  = 1,2,3)  allows  us  to  write  \j0  as 

V0  = (v0)jfj  (20) 


in  which  (y 0)j  = V0  • f )■ 
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Suppose  that  v^Dy^  and  D2y0  are  linearly  dependent  for  all  D and 
V$,  there  are  a,  (3  and  7 which  are  not  all  zero,  such  that 

asjO  + PD  v#  + 7D2  v#  = 0.  (21) 

Substituting  (19)  and  (20)  into  (21)  yields 

3 

5^(a  + /9/ifc  + 7/i|)(v5)fcffc  = 0,  (22) 

k= 1 

which  implies,  as  fk  ( k — 1,2,3)  are  linearly  independent, 

(a  + pfik  + ^l)(v0)k=0,  (k  = 1,2,3).  (23) 

For  arbitrary  ^/O,  (v^)fe  need  not  be  zero,  so 

+ j9/ifc  + 7/4  = 0,  (k  = 1,2,3)  (24) 

that  requires  that  a = (3  = 7 = 0 for  distinct  fxk,  contrary  to  the  hypothesis. 
Theorem  4 has,  thus,  been  proved. 

Applying  Theorem  4 to  the  SGS  flux  vector,  we  have 

f (0,  o.p,  V0,  d)  = 0O(0,  o.p,  v0,  d)  + 01  (0,  o.p,  v0»  d)d 
+02(0,O.P,V0,D)D2V0, 

(25) 

and 

f (0, 0.P.,  - V0,  D)  = -0O(0, OP,  - V0,  D)  V0 

-01  (0,  O.P,  - V0,  D)D V0  - 02 (0, OP,  - v0,  D)D2  V0.  ^ ' 

Substituting  (25)  and  (26)  into  (18)  leads  to 

[0o (0, 0.P.,  D)  - 0o (0,  O.P,  - v0,  D)] y0 


d" [01.(0,  O.P,  V0,  D)  - 0i(0,  O.P,  - V0,  D)]DV0  (27) 

+[02(0,  O.P,  V0,  D)  - 02(0,  O.p,  - V0,  D)]D2  V0  = 0 

which  implies,  since  and  D2v0  are  linearly  independent, 

0i(0,O.P,v0,D)  = 0i(0,O.P,-v0,D),  (*  = 0,1,2).  (28) 

To  satisfy  this  requirement,  take 

0i(0,  O.P,  v0>  D)  = 0i(0,  O.P,  y0  ® V0>  D),  (*  = 0,1,2).  (29) 

Then  (25)  and  (29)  result  in 

f (0,  O.P,  QV0,  QDQr)  = Q(0O  V0  + 0iD  V0  + 02D2  V0)  (30) 
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in  which, 

A = A(0, 0.P.,  Qv#  ® Qv#,  QDQt), 

and 

Qf(0, O.P.,  s/d,  D)  = Q(</>0  V#  + ^iDv0  + ?/>2D2v0)-  (31) 

By  making  use  of  (30)  and  (31),  (17)  yields 

W>o  - + {A  - + {i>2  - A)P> 2V0  = 0 (32) 

that  implies,  by  Theorem  4, 

^(0,O.P,QV0®QV0,QDQT)  = ^(0,O.P,V0®V^D)  VQ.  (33) 
Theorem  5.  Suppose 

ip(6, 0.P.,  Qb  ® Qb,  QBQt)  = ip(9, 0.P.,  b®b,  B),  Vb  and  B, 

then 

ip(6, 0.P.,  a ® a,  A)  = ip(6,O.P.,  b ®b,  B) 

whenever  Jk(a,  A)  = Jk{b,B)  (k  = 1,2,  ••■,6).  Here 

Ji(a,  A)  = trA,  J2(a,A)  = ^[(frA)2  - fr(A2)],  a,  A)  = detA, 
J4(a,  A)  = a • Aa,  J5(a,  A)  = a • A2a,  Je(a,  A)=|a|, 

a and  b are  two  arbitrary  vectors,  A and  B are  two  arbitrary  symmetric 
tensors. 

Proof  Since  Jk{ a,  A)  = Jk{ b,  B)  ( k = 1,2, 3),  tensors  A and  B have  same 
eigenvalues.  Let  pk  be  their  eigenvalues,  A and  B may  be  written  as, 

3 3 

A = 'y  ^ Pk&k  ® ®fc?  B = y ) Pk^k  ® ffc 

fc= 1 fc=l 

where  ek  and  ffc  (k  = 1,2,3)  are  eigenvectors  of  A and  B,  respectively. 
Define 

Q — €fc  G1  f fc 

that  is  a rotation  tensor,  and 

e*  = Qfp  A = QBQt,  A2  = QB2Qt,  (i  = 1,2,3).  (34) 

By  applying  Jk( a,  A)  = Jk( b,B)  ( k - 4,5,6),  we  have 

3 3 3 3 

X^(b  • ffc)2  = X)(Qra  • ffc)2,  yr  Pk{ b • ffc)2  = Y Mfe(QTa  • ffe)2, 

k=  1 k= 1 k= 1 k=  1 

Y dl(h  • ffc)2  = Y dl(QTa  ■ ffc)2- 

fc=i  fc=i 


(35) 
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This  implies,  for  the  distinct  p,k  (k  = 1, 2, 3), 

(QTa  b)  • f*;  = 0,  (k  = 1,2,3).  (36) 

Note  that  fj,  ( k = 1, 2, 3)  are  linearly  independent,  then  a = ±Qb,  a® 
a = Qb  0 Qb.  By  hypothesis, 

0.P. , b 0 b,  B)  = 0(0, 0.P.,  Qb  0 Qb,  QBQT)  = 0(0, 0.P.,  a 0 a,  A) 
in  which  A = QBQr  [ (34)]  and  a 0 a = Qb  0 Qb  are  used.  Therefore, 
#JO.P,b«b,B)=#,O.P,Jk(b,B)],  (k  = 1,2,  •••,6)  (37) 

if 


0(0, 0.P.,  Qb  0 Qb,  QBQr)  = 0(0,  O.P.,  b 0 b,  B),  Vb  and  B.  (38) 

The  converse  is  also  true  since  Jk( Qb,  QBQr)  = Jfc(b,  B)  (k  = 1, 2,  • • • , 6). 

Theorem  6.  The  necessary  condition  for  the  constitutive  process  (1)  to 
satisfy  requirements  (2)  and  (3)  is 

q = f (6, 0.P.,  V0,  D)  = (</)0 1 + 0iD  + 02D2) 


where 


(pi  = 4>i[0,O.P.,  J*;(v0,I>)],  (*  = 0,1,2;  k = 1,2,  ■ • • ,6). 

Proof  Applying  Theorem  5 to  (33)  yields 

0i(0, 0.P.,  v0  <8>  V*.  D)  = 0*[0, 0.P.,  Jfc(v^,  D)].  (39) 

This,  with  (25)  and  (29),  leads  to 

q = f(0, 0.P.,  vfl,  D)  = (0O1  + 0iD  + 02D2)  (40) 


where 


4>i  = ^[9, 0.P. , Jfc(V0,  D)],  (<  = 0, 1, 2;  k = 1, 2,  • • • , 6). 

If  the  three  eigenvalues  of  D are  not  distinct,  we  can  still  obtain  (40)  with 
01  = 02  = 0 (for  the  case  of  three  identical  eigenvalues)  or  02  = 0 (for  the 
case  of  two  identical  eigenvalues)  by  the  similar  method.  Therefore,  (40)  is 
valid  for  all  cases. 
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2.2.  SUFFICIENCY  OF  (40)  FOR  REQUIREMENTS  (2)  AND  (3) 

Suppose  (40)  holds,  then 

f (0*,  o.p.*,  (V0*)*,v*,L*)  = {M0\o.ps,  Jfc((v*T,D*))i 
+01(r,O.PA,Jfc((v^)*,DQ)D*  + </»2(r,O.PA,Jfc((Vr)*,D*))D*2} 
MT  - {0o(<9,  O.p,  J,(Q  v QDQt))Q1Qt  + MO,  O.P. , Jfc(Q  V 
QDQt))QDQt  + 02(0,  O.P,  Jk{ Q V 0,  QDQT))QD2Qr}Q  V e 
= Q{M0 , O.p,  4(Q  V 0,  QDQt))1  + 0!  (0,  O.P,  Jfc(Q  V 0,  QDQt)) 

D + MO,  O.P,  J*(Q  V 0,  QDQt))D2}QtQ  v 0 
= Qf (0,  O.P,  V0,  V,  L)  = q*,  (A  = 1, 2,  ■ • • , 6) 

(41) 

iii  which  Eq.(5),  D*  = QDQ7  and  Jk( Q v 0,  QDQT)  = Jfc(y0,D)  ( k = 

1. 2,  • • ■ , 6)  are  used. 

Also,  if  (40)  holds, 


f (0,  O.P,  - v 0,  v,  L)  = {MO,  O.p.,  Jk{-  V 9,  D))l 

+0i ( 9 , O.P,  J*(-  V 9 , D))D  + 02(0,  O.P,  Jfc(-  V 0, D))D2}(—  V 0) 

= — {0O(0,  OP,  Jfc(V0,  D))l  + 0i (0,  OP,  Jk(s/0 , D))D 

+02(0,  O.P,  Jfe(V0,  D))D2}(V0)  = -f(0,  O.P,  V0,  v,  L),  (fc  = 1, 2,  • • • , 6) 

(42) 

in  which  Jk(—  y0,D)  = J^(y0,D)  ( k = 1,2, •••,6)  are  used.  Equations 
(41)  and  (42)  establish  the  sufficiency  of  (40)  for  (2)  and  (3). 

2.3.  PROPERTIES  OF  0/  (/  = 0,1,2)  AND  BOTH  NECESSARY  AND 
SUFFICIENT  CONDITIONS  FOR  INEQUALITY  (4) 

Both  necessity  and  sufficiency  of  (40)  for  Eqs.(2)  and  (3)  are  established 
in  §2.1  and  §2.2.  Here  we  analyze  some  fundamental  properties  of  0j  (i  = 
0, 1,  2),  and  develop  both  necessary  and  sufficient  conditions  for  inequality 

(4). 

Rewrite  Eq.(40)  as 

q = -Ky0,  (43) 

with 

K = (0o  1 + 0i  D + 02D2)  = <p0l  + (piD  + <p2  D2.  (44) 

Here  1 is  a unit  (identity)  tensor,  and 

<Fi  = "0i  (*  = 0,1,2).  (45) 


K has  two  fundamental  properties:  (1)  it  is  a real-valued  tensor  on  the 
ground  of  practical  transport  processes,  and  (2)  it  is  a symmetric  tensor 
due  to  the  symmetry  of  velocity  strain  tensor  D. 
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Let  Xj  and  fj  {j  — l, 2, 3)  be  the  three  eigenvalues  of  D and  K,  respec- 
tively. Since  K is  related  to  D through  Eq.(44), 

fj  = <fo  + <pi\j  + ( j = 1, 2, 3).  (46) 

Because  D and  K are  real-valued  symmetric  tensors,  A j and  fj  ( j = 1, 2, 3) 
must  be  real- valued,  i.e., 


fj  = fj,  (j  = 1,2,3), 

(47) 

Xj  — Xj , ( j = 1,2,3). 

(48) 

By  Eq.(46), 

7 j = Vo  + Xj  + ^PyXj , (j  = 1,2,3). 

(49) 

By  making  use  of  Eqs.(47)  and  (48),  Eq.(49)  leads  to 

fj  = ipo  + ViXj  + <p2X2j,  (j  = 1,2,3). 

(50) 

This,  with  Eq.(46),  yields 

(<A)  ~ ¥>o)  + (^l  - Vi)Xj  + {tfi2  ~ <P2)tf  = °>  U = !,  2, 3), 

VA?  g R 

(51) 

which  indicates  that 

<Pi  = Vi  (*  = 0,1,2). 

(52) 

Therefore,  </?,  (i  = 0, 1, 2)  must  be  real-valued. 
Substituting  Eq.  (43)  into  inequality  (4)  yields 


V#  • K v#  > 0,  Vv#,  (53) 

which  implies  that  K is  positive  semi-definite.  Note  also  that  K is,  in  prac- 
tice, an  invertible  tensor,  it  must  be  positive  definite.  The  same  conclusion 
may  be  obtained  by  noting  that  the  equal  sign  in  (4)  is  only  for  reversible 
processes  and  transport  processes  are  irreversible. 

The  necessary  and  sufficient  condition  for  a symmetric  tensor  to  be 
positive  definite  is  that  all  of  its  eigenvalues  are  positive  definite.  Both 
necessary  and  sufficient  condition  for  inequality  (4)  is,  thus, 

¥?o  + ViXj  + <p2Xj  > 0,  ( j = 1, 2, 3),  VA j 6 R.  (54) 

Two  necessary  conditions  of  (54)  can  be  easily  obtained  by  considering 
cases  of  A j — 0 and  | Xj  |— > oo,  respectively,  as 


(fo  > 0, 


(55) 


PHYSICS-PRESERVING  TURBULENT  CLOSURE  MODELS 


273 


<P2  > 0.  (56) 

Dividing  (54)  by  tpQ.  we  can  rearrange  (54)  into  an  alternative  form 

(1  + ^)2  + (^  - ^\)\)  > 0 VA j e R.  (57) 

2<p0  Vo  W 3 J 

This  yields  another  necessary  condition,  by  setting  A j = — 2<po/Vi> 

ipf  - 4ip0(p2  < 0.  (58) 

Conversely,  it  is  easy  to  show  that  (55),  (56)  and  (58)  are  also  the  sufficient 
conditions  of  (54). 

The  detailed  expressions  of  (po,<Pi  and  ip2  are  material-dependent  and 
need  to  be  determined  through  experiments.  Once  they  are  determined, 
Eq.  (43)  can  serve  as  the  SGS  flux  model  that  is  properly  invariant  and 
satisfies  the  second  law  of  thermodynamics. 

3.  Concluding  Remarks 

For  a class  of  turbulence  flows  for  which  the  SGS  flux  vector  can  be  de- 
scribed by  Eq.(l),  both  necessary  and  sufficient  conditions  are  derived  in 
a systematic,  rigorous  way  for  the  invariance,  the  PMFI  and  the  second 
law  of  thermodynamics.  This  leads  to  a general  model  (43)  with  three  real- 
valued functions  <pi  ( i = 0, 1,2)  satisfying  (55)  (56)  and  (58).  Any  specific 
model  satisfying  (43),  (55),  (56)  and  (58)  is  properly  invariant  and  satisfies 
the  second  law  of  thermodynamics  while  the  model  violating  these  condi- 
tions is  not.  The  work  is  believed  to  be  important  both  for  developing  the 
specific,  physics-preserving  models  and  for  clarifying  some  confusion  in  the 
literature  by  noting  that  the  previous  works  employ  an  intuitive  approach 
with  the  focus  on  only  obtaining  the  sufficient  condition. 
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